WKB ANALYSIS FOR NONLINEAR SCHRODINGER 
EQUATIONS WITH POTENTIAL 



REMI CARLES 

Abstract. We justify the WKB analysis for the semiclassical nonlinear 
Schrodinger equation with a subquadratic potential. This concerns sub- 
critical, critical, and supercritical cases as far as the geometrical optics 
method is concerned. In the supercritical case, this extends a previous 
result by E. Grenier; we also have to restrict to nonlinearities which are 
defocusing and cubic at the origin, but besides subquadratic potentials, 
we consider initial phases which may be unbounded. For this, we con- 
struct solutions for some compressible Euler equations with unbounded 
source term and unbounded initial velocity. 



1. Introduction 
Consider the initial value problem, for x 6 M. n and k > 0: 

(1.1) ied t u £ + —Au £ = V(t, x)u £ + e K f (\u £ \ 2 ) u £ 

(1.2) u £ t=0 = flg(«)^«/«. 

The aim of WKB methods is to describe u £ in the limit e — > 0, when <f>o 
does not depend on e, and Oq has an asymptotic expansion of the form: 

(1.3) a 6 o( x ) ~ ao{x) + ea\{x) + e 2 a 2 (x) + . . . 

The parameter k > describes the strength of a coupling constant, which 
makes nonlinear effects more or less important in the limit e —* 0; the larger 
the k, the weaker the nonlinear interactions. In this paper, we describe the 
asymptotic behavior of u £ at leading order, when the potential V and the 
initial phase (fio are smooth, and subquadratic in the space variable. 

Such equations as 1)1. 1|) appear in physics: see e.g. [HI] for a general 
overview. For instance, they are used to model Bose-Einstein condensation 
when V is an harmonic potential (isotropic or anisotropic) and the non- 
linearity is cubic or quintic (see e.g. ^1 E3 EDI)- In most of this paper, 
the initial data we consider are in Sobolev spaces H s . We refer to |Hj for 
numerics on the semi-classical limit of (|1.1|) . 

We shall not recall the results concerning the Cauchy problem for ()1.1|) - 
(|1.2|) . and refer to jlUj for an overview on the semilinear Schrodinger equa- 
tion. 
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In the one-dimensional case n = 1, the cubic nonlinear Schrodinger equa- 
tion is completely integrable, in the absence of an external potential, or when 
V is a quadratic polynomial (HI p. 375]). Several tough papers analyze the 
semi-classical limit in the case V = 0, for k = 0: see e.g. [25*1 l2"""l "3l"l "*"""] . 
We shall not use this approach, but rather work in the spirit of |22j . 

An interesting feature of (jl.lj) is that one does not expect the creation 
of harmonics, provided that only one phase is present initially, like in (|1.2j) . 
The WKB methods consist in seeking an approximate solution to of 
the form: 

(1.4) u £ (t, x) ~ (ao(t, x) + eai(t, x) + e 2 a 2 (t, x) + . . .) e**^*)/ 6 . 

One must not expect this approach to be valid when caustics are formed: 
near a caustic, all the terms <£, ao, ai, . . . become singular. Past the caustic, 
several phases are necessary in general to describe the asymptotic behavior 
of the solution (see e.g. |17j for a general theory in the linear case). In this 
paper, we restrict our attention to times preceding this break-up. 

For such an expansion to be available with profiles a.j independent of e, 
it is reasonable to assume that k is an integer. However, we do not assume 
that k is an integer: we study the asymptotic behavior of u e at leading order 
(strong limits in L 2 f] L°° for instance), including cases where other powers 
of e would come into play. 

We distinguish two families of assumptions: "geometrical" assumptions, 
on the potential V and the initial phase cfto, and "analytical" assumptions, 
on / and the initial amplitude Oq. In all the cases, we shall not try to seek 
the optimal regularity; we focus our interest on the limit e — ► 0. 

Assumption 1 (Geometrical). We assume that the potential and the initial 
phase are smooth and subquadratic: 

• V G C°°(M t x K£), and d*V G L™ c (R t ; L°° (R™)) as soon as \a\ > 2. 

• (Ao G C°°(R n ), and Sg^ G L°°(M n ) as soon as \a\ > 2. 

The assumption of V being subquadratic is classical in other contexts; 
for instance, locally in time, the dispersion for e~ l ~^~ z A+v ^ is the same as 
without potential (see [TU1 [2*0] ). 



hence yielding the same local Strichartz estimates as in the free case. Global 
in time Strichartz estimates must not be expected in general, as shown by 
the example of the harmonic oscillator, which has eigenvalues. For positive 
superquadratic potentials, the smoothing effects and Strichartz estimates 
are different (see (HI1I38 ). This is related to the properties of the Hamilton 
flow, which also imply: 



Lemma 1.1. Under Assumption^ there exist T > and a unique solution 
4ik G C°°([0,T] x R n ) to: 

(1.5) d t (t> cik + ^iV^cikl 2 + V(t,x) = ; 4>eik\t=o = <t>o ■ 

This solution is subquadratic: <9"</> e ;k G L°°([0,T] x M n ) as soon as \a\ > 2. 
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This result is proved in Section |21 where other remarks on Assumption ^ 
are made. 

Assumption 2 (Analytical). We assume that the nonlinearity is smooth, 
and that the initial amplitude converges in Sobolev spaces: 

• / G C°°(M;M). 

• There exists ao G H°° := Pi s >oH s (W n ) , such that aff converges to ao 
in H s for any s, as e — ► 0. 

Remark. Some of the results we shall prove remain valid when / is complex- 
valued. In that case, the conservation of mass associated to the Schrodinger 
equation, [|u e (t)[|£2 = ||«oI|l 2 ) no longer holds. On the other hand, when 
< k < 1, this assumption is necessary in our approach, and we even assume 
/' > 0. 

1.1. Subcritical and critical cases: k > 1. When the initial data is of 
the form (jl..3j) . the usual approach consists in plugging a formal expansion 
of the form 1)1.4)1 into 1)1.1)) . Ordering the terms in powers of e, and canceling 
the cascade of equations thus obtained yields <3?, ao, ai, ... 

Assume in this section that k > 1, and apply the above procedure. To 
cancel the term of order O(e ), we find that $ must solve 1)1.5)) : = ipd^. 
Canceling the term of order 0(e 1 ), we get: 



1 [o 

<9ta + Vc/>eik • Va + -a A^ cik = <^ 

2 -*/(|ao| 



if k > 1, 
if as = 1. 



We see that the value k = 1 is critical as far as nonlinear effects are con- 
cerned: if k > 1, no nonlinear effect is expected at leading order, since 
formally, u e ~ aoe l ^ cik//e , and </> e ik and ao do not depend on the nonlinearity 
/. If k = 1, then ao solves a nonlinear equation involving /. 

We will see in Section El that ao solves a transport equation that turns out 
to be a ordinary differential equation along the rays of geometrical optics, 
as is usual in the hyperbolic case (see e.g. [SB)- More typical of Schrodinger 
equation is the fact that this ordinary differential equation can be solved 
explicitly: the nonlinear effect is measured by a nonlinear phase shift (see 
the example of j2jj for a similar result in the hyperbolic setting). We prove 
the following result in Section |2J 

Proposition 1.2. Let Assumptions^ and\^ be satisfied. Let k > 1. Then 
for any e g]0, 1], $£$-$12$) has a unique solution u £ £ C°°([0,T] x R n ) n 
C([0, T];H S ) for any s > n/2 (T is given by Lemm,a, Moreover, 
there exist a,G G C°°([0, T] x 1"), independent of e G]0, 1], where a G 
C([0,T];L 2 n L°°), and G is real-valued with G G C([0, T];L°°), such that: 



u £ -ae ieK lG e^ eik / £ 



L°°([0,r];L 2 nL°°) 



as e 



0. 



The profile a solves the initial value problem: 

1 



(1.6) d t a + V0 cik • Va + -aa<p eik = u ; a\ t=0 = a , 

and G depends nonlinearly on a through f. In particular, if k > 1, then 



L<x>([0,T];L 2 nL°°) 



as e — > 0, 
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and no nonlinear effect is present in the leading order behavior of u £ . If 
k = 1, nonlinear effects are present at leading order, measured by G. 

The dependence of G upon a and / is made more explicit in Section [3 
in terms of the Hamilton flow determining <^ e ik (see (|3.4|l ). Note that in the 
above result, we do not assume that k is an integer. 

1.2. Supercritical case: k = 0. It follows from the above analysis that 
the case < k < 1 is supercritical. We restrict our attention to the case 
k = 0. We present an analysis of the range < k < 1 in Section H3 Plugging 
an asymptotic expansion of the form (|1.4|) into yields a shifted cascade 
of equations: 

O(e ): a^ + i|V$| 2 + y + /(|a | 2 ) =0, 

O {e 1 ) : a t a + V$ • Va + ^a A$ = 2if (|a | 2 ) Re (a al) . 

Two comments are in order. First, we see that there is a strong coupling 
between the phase and the main amplitude: ao is present in the equation for 
Second, the above system is not closed: <3? is determined in function of 
ao, and ao is determined in function of ai. Even if we pursued the cascade 
of equations, this phenomenon would remain: no matter how many terms 
are computed, the system is never closed (see (25 )• This is a typical feature 
of supercritical cases in nonlinear geometrical optics (see j!21 113) ) . 

In the case when V = and cfto £ H s , this problem was resolved by 
E. Grenier [22], by modifying the usual WKB methods; this approach is 
recalled in Section 0J Note that even though ai is not determined by the 
above system, the pair (p,v) := (|ao| , V3>) solves a compressible Euler 
equation: 

dtv + v- Vv + VV + Vf(p) =0 ; v\ = V0 O 
(1.7) , 

d t p + V • (pv) = ; p\ t=Q = \a \ . 

Using techniques introduced in the study of quasilinear hyperbolic equa- 
tions, E. Grenier justified a WKB expansion for nonlinearities which are 
defocusing, and cubic at the origin (/' > 0). We shall not change this as- 
sumption, but show how to treat the case of a subquadratic potential with 
a subquadratic initial phase. Note that even the construction of solution 
to (|1.7|) under Assumption ^ is not standard: the source term W may be 
unbounded, as well as the initial velocity V0o- 

Assumption 3. In addition to Assumption^ we assume: 

• f>0. 

• There exists «o> a i £ H°° , with xa$,xai G H°° , such that: 
||«o ~ a o ~ £a i\\H s + ll xa o — xa o ~ sxai\\ Hs = o(e), Vs > 0. 

We can then describe the asymptotic behavior of the solution to (jl.lj) - 
1)1.2(1 for small times: 

Theorem 1.3. Let Assumptions^ and\^ be satisfied. Let k = 0. There 
exists T* > independent of e e]0, 1] and a unique solution u e S C°°([0, T#] x 
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R n ) n C([0,T*};H s ) for any s > n/2 to ifTTJl -lfT^ ji . Moreover, there exist 
a,(p £ C([0, T*]; H s ) for every s > 0, such that: 

(1.8) limsup u £ - ae i( * + ^ ik)/£ = 0(t) as t -» 0. 



£^0 



L 2 nL°° 



Here, a and (j) are nonlinear functions o/0 e ik andao, given by 1)5.9 J) . Finally, 
there exists e C([0, T*]; .£P) /or every s > 0, real-valued, such that: 



(1.9) limsup sup 

e^O 0<t<T* 



ae i</> (1) e «(<A+</'oik)/£ 



= 0. 

L 2 nL°° 



JTie phase shift (f)^ 1 ' is a nonlinear function o/^ e ik,«o a i- 

This result can be understood as follows. At leading order, the amplitude 
of u £ is given by ae l ^ (1) , which can be approximated for small times by a, 
because 0^L_ O = 0. The rapid oscillations are described by the phase 
<t> + '/'oik- The function <\> is constructed as a perturbation of <^> e ik) but must 
not be considered as negligible: its .£P-norms are not small in general, see 
1)5.9)1 (at time t = for instance). As a consequence of our analysis, the pair 



(p,v) = (|a| 2 ,V(0 + eik )) 
solves the system 1)1.7)1 . 



ae^ (1) 



2 

, V(0 + </>eik! 



Remark 1.4. With this result, we could deduce instability phenomena for 
1)1.1)1 - 1)1.2)1 in the same fashion as in 0. Note that because of Assumption 
it seems that the approaches of [£l E] cannot be adapted to the present 
case: the Laplacian can never be neglected, and apparently, WKB approach 
is really needed. 

The rest of this paper is organized as follows. In Section [21 we prove 
Lemma ll. II In Sectional we prove Proposition 11.21 and explain how G is 
obtained. We recall the approach of [22] in Section [3J and show how to 
adapt it to prove Theorem 11.31 in Sectional We present an analysis for the 
case < k < 1 in Section [J3 

2. Global in space Hamilton-Jacobi theory 
In this section, we prove Lemma ll.ll Consider the classical Hamiltonian: 

H(t, x, r, £) = r + ^|£| 2 + V(t, x), (t, i,T,Qel + xR"xRx R n . 

It is smooth by Assumption^ Therefore, it is classical (see e.g. jlHj ) that in 
the neighborhood of each point x £ W 1 , one can construct a smooth solution 
to the eikonal equation 1)1.5)) . on some time interval [—t(x),t(x)}, for some 
t(x) > depending on x. The fact that in Lemma ll.l) we can find some 
T > uniform in x £ M. n is due to the fact that the potential and the initial 
phase are subquadratic. 

Recall that if there exist some constants a, b > such that a potential V 
satisfies V(x) > — a\x\ 2 — b, then — A + V is essentially self-adjoint on Cq ^™) 
(see [SH p. 199]). If —V has superquadratic growth, then it is not possible to 
define e - i '(- A+v ) (see QUI Chap. 13, Sect. 6, Cor. 22] for the case V(x) = -x 4 
in space dimension one). This is due to the fact that classical trajectories 
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can reach an infinite speed. We will see below that if the initial phase 4>q 
is superquadratic, then focusing at the origin may occur "instantly" (see 
Example ^) . 

To construct the solution of (|1.5j) . introduce the flow associated to H: let 
x(t,y) and £(t,y) solve 

J d t x{t,y) =£(t,y) ; x(0,y)=y, 
\d t ((t,y) = -V x V(t,x(t,y)) ; £(0, y) = Vfo(y). 

Recall the result (valid under weaker conditions than Assumption : 

Theorem 2.1 (|16j. Th. A. 3. 2). Suppose that Assumption^ is satisfied. Let 
t G [0, T] and 9q an open set ofM. n . Denote 

9 t :={x(t,y)eR n ,y£9 } ; 9 := {(t, x) G [0, T] x R n , x € 0*}. 

Suppose that for t G [0, T], t/ie mapping 

9 B y^ x(t,y) G t 

is bijective, and denote by y(t, x) its inverse. Assume also that 

V x y G 

T/ien i/iere exists a unique function 9 3 (t,x) t— > iji>eik(^ ^) £ K solves 
1)1.5(1 . and satisfies V^ e ik G Lj^ c (0). Moreover, 

(2.2) V^eik(t,ar) = e(i,y(i^))- 

Proposition 2.2 ([33], Th. 1.22 and QS|, Prop. A.7.1). Suppose that the 
function M. n 3 y \— * x{y) G M n satisfies: 

| det V y x| > C > and |fl£x| < C, |a| = 1, 2. 

T/ien x is bijective. 

Proof of Lemma M.lX Lemma 1 1 . 1 1 follows from the above two results. From 
Assumption^ we know that we can solve (|2.1jl locally in time in the neigh- 
borhood of any y G M n . Differentiate (|2.1|) with respect to y: 

j d t d y x{t, y) = d y i {t, y) ; d y x(0, y) = Id, 

[d t d y at,y) = -v 2 x v(t,x(t,y))d y x(t,y) ■ d y ao,y) = v 2 My)- 

Integrating 1)2.3(1 in time, we infer from Assumption U that for any T > 0, 
there exists C T such that for (t,y) G [0,T] x R n : 

\d y x(t,y)\ + \d y C(t,y)\<C T + C T f (\d y x( S ,y)\ + \dyt(s,y)\)ds. 

Jo 

Gronwall lemma yields: 

(2-4) \\d y x(t)\\ LX + \\d y at)\\L^<C'(T). 

y y 

Similarly, 

(2.5) || W)|| x?0 + \\d^(t)\\ L ^ < C(a,T), Va G TC, \a\ > 1. 
Integrating the first line of 1)2.3(1 in time, we have: 

det V y x(t, y) = det ^Id + ^ V<,£ (s, y) ds^j . 
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We infer from Q2.4|) that for t G [0, T], provided that T > is sufficiently 
small, we can find Cq > such that: 

(2.6) \detV y x(t,y)\ > C , V(i,y) G [0,T] x K n . 

Applying Proposition 12/21 we deduce that we can invert y i— > x(t, y) for 
t€ [0,T]. 

To apply Theorem 12.11 with 6q = 8 = 9t = K n , we must check that 
V x y G L^ c (lR n ). Differentiate the relation 

x(t,y(t,x)) = x 

with respect to x: 

V x y(t, x)V y x (t, y(t, x)) = Id. 
Therefore, V x y(t,x) = V y x (t, y(t, as matrices, and 

(2.7) V x y(t, x) = dety l x{s -jy adj (V y x (t, y(i, a;))) , 

where adj (V^x) denotes the adjugate of V y x. We infer from (|2,4j) and 
(t2~fij) that V^y G L°°(]R n ) for t G [0,T]. Therefore, Theorem [2j] yields 
a smooth solution </> e ik 

to (|1.5j) . local in time and global in space: </> e ik £ 

C°°([0,T] x R n ). 

The fact that </> e ik is sub quadratic as stated in Lemma 11.11 then stems 
from D, (|231), C02> and <22J. □ 



We now give some two examples showing that Assumption ^ is essentially 
sharp to solve IJ1.5j) globally in space, at least when no assumption is made 
on the sign of V nor on the geometry of V0o- We already recalled that if 
— V has a superquadratic growth, then —A + V is not essentially self-adjoint 
on C^°(IR n ), so we shall rather study the dependence of ^> e ik on the initial 
phase 4>q. 

Example 1. Assume that V = and 

Mx) = ~ {2 + \ 5)T {\x\ 2 + l) 1+S , T>0, 5^-1. 

Then Assumption^ is satisfied if and only if 8 < 0. When 5 = 0, then (|1.5() 
is solved explicitly: 

\x\ 2 1 

<Peik(t,x) 



2(t-T) 2T' 

This shows that we can solve (|1.5(l globally in space, but only locally in time: 
as t — > T, a caustic reduced to a single point (the origin) is formed. Note 
that with T < 0, 1)1.5(1 can be solved globally in time for positive times. 
When 5 > 0, then integrating (|2.1j) yields: 

%{t, y) = y + J £0> y)ds = y + J £(0, y)ds = y - ^ (|y| 2 + l)* 5 y 



\y\ 2 + i) s ) 
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For R > 0, we see that the rays starting from the ball {\y\ = R} meet at 
the origin at time 

T 

Tc{R) = (R? + iy 

Since R is arbitrary, this shows that several rays can meet arbitrarily fast, 
thus showing that Theorem 12. II cannot be applied uniformly in space. 

Example 2. When V(t, x) = ^ Y^j=\ ^f^f * s an harmonic potential (u>j / 0), 
and (fro = 0, we have: 

n 

(frit, x) = - y x i tan(u^i). 
i=i 

This also shows that we can solve 1)1. 5JI globally in space, but locally in time 
only. Note that if we replace formally ujj by iuj, then V is turned into —V, 
and the trigonometric functions become hyperbolic functions: we can then 
solve p.5j) globally in space and time. 

Instead of invoking Theorem 12.11 and Proposition 12.21 one might try to 
differentiate ()1.5|) in order to prove that (fr^ is subquadratic, in the same 
fashion as in |1J|H]. For 1 < j, k < n, differentiate (|1.5|) with respect to Xj 
and Xk- 

n 

d t dj k 4> cik + V x 0eik • V x (dj k (fr cik ) + ^2 dji4>dkd? k (freik + dj k V(t, x) = ; 

i=i 

<9ffc0eik|t=O = 9j k <fro ■ 

We see that we obtain a system of the form 

D t y = Q(y)+R ; y\ t=Q = y(0), 

where y stands for the family (<9| fc <?i>eik)i<.7,fc<n> Q is quadratic, and R and 
y(0) are bounded. The operator Dt is a well-defined transport operator 
provided that the characteristics given by: 

d t x(t,y) = V x (fr cik (t,x(t,y)) ; x(0,y)=y, 

define a global diffeomorphism. Proving this amounts to using Proposi- 
tion 12.21 for the rather general initial data we consider. So it seems that 
this approach does not allow to shorten the proof of Lemma 11.11 

3. SUBCRITICAL AND CRITICAL CASES 
To establish Proposition 11.21 define 

a £ (t,x) :=u £ (t,x)e- i ^ {t ' xy£ . 
Then u £ solves (|1.1|) - (|1,2|) if and only if a £ solves: 

d t a £ + V(fr cik ■ Va £ + \a £ A(fr e]k = i £ -Aa £ - ie"" 1 / (|a e | 2 ) a £ , 
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Proposition 3.1. Let Assumptions^ and\^be satisfied. Let k > 1. For any 

e e]0,l], (E3) fcaa a umgue sotoion a £ £ C°°([0,T] x M n ) n C([0, T}; H s ) 
for any s > n/2. Moreover, a £ is bounded in L°°([0,T];H s ) uniformly in 
£ G]0, 1], for any s > 0. 

Proof. Using a mollification procedure, we see that it is enough to establish 
energy estimates for (f3. 1|) in H s , for any s > 0. Let s > n/2, and a E N n , 
with s = |a|. Applying to (|3.1j) . we find: 

(3.2) d t d*a £ + V0 eik • V0£a e = i|Ad> £ - iz K ~ x d* (f (\a £ \ 2 ) a £ ) + R e a , 
where 



^ = [V4ik"V,^] a £ --«9«(a £ 



eikj 



Take the inner product of (|3.2j) with d a a £ , and consider the real part: the 
first term of the right hand side of (|3.2j) vanishes, and we have: 



2dt u x 



L-' 



+ Re / d%a £ V^ eik • V<9> £ ^e*" 1 \\f {\a £ \ 2 ) a £ \\ Hs \\a £ \\ HS 



+ \\R 



Notice that we have 

Re / d«a7V<p eik -Vd°a £ 



a\\L 2 



V</»eik-V|a> £ | 



d°a £ \ 2 A^ eik 



< C \\d 



£||2 



since A^ e jk € L°°([0,T] x 1") from Lemma ll.ll Moser's inequality yields: 

||/(KlV||„ s <c(|K| 

Summing over a such that \a\ = s, we infer: 
d 



dt 



\a £ \\ H s < C(||a e || LOO ) \\a £ \\ HS + \\R 



a\\H s ■ 



Note that the above locally bounded map C(-) is independent of e if and 
only if k > 1. To apply Gronwall lemma, we need to estimate the last term: 
we use the fact that the derivatives of order at least two of </> e ik are bounded, 
from Lemma ll.l| to have: 

We can then conclude by a continuity argument and Gronwall lemma: 

IKIU°°([0,T];.fP) < C (s, HoqII^.) • 

For the mollification procedure, this yields boundedness in the high norm. 
Contraction in the small norm then follows easily with classical arguments 
(see e.g. OEH]), completing the proof of Proposition 13.11 □ 

Corollary 3.2. Let Assumptions^ and\^ be satisfied. Let k > 1. Then 



\L°°([0,T};H s ) 



as e -> 0, Vs > 0, 



where a £ solves: 



(3.3) d t a £ + V^ cik • Va £ + -a £ 



'eik 



-IE 



K-l 



f(\a 



; a \t=o - a o- 
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Proof. The proof of Proposition EHl shows that a £ £ C°°([0,T] x R n ), and 
that a £ is bounded in L°°([0,T]; H s ) uniformly in e £]0, 1], for any s > 0. 
Let w £ = a £ — a £ : it solves 

<9 iU ; e + • Vw £ + ^ £ A0 cik = i e -Aa £ - ie^ 1 (F(a £ ) - F(a £ )) , 

wft=o = a o - a o, 

where we have denoted -F(.z) = /(|z| 2 )z. Proceeding as in the proof of 
Proposition 13.11 we have the following energy estimate: 

4lKH.ffs < C\\w £ \\ Hs +e\\Aa £ \\ Hs + \\F(a £ ) - F(a £ )\\ Hs . 

Since H s is an algebra and F is C , the Fundamental Theorem of Calculus 
yields: 

||F(a £ ) - F(a £ )\\ Hs < C (\\a £ \\ Hs , \\a £ \\ Hs ) \\w £ \\ H s. 
Now since a £ and a £ are bounded in L°°([0, T]; iF+ 2 ) uniformly in e £]0, 1], 
we have an estimate of the form: 

4lKllff s < C(s) ||w E ||^ + C(s)e. 

We conclude by Gronwall lemma, since 1 1 a,g — a-o II J/ s ~^ from Assumption |2 

□ 

Remark 3.3. If we assume moreover that like in (|1.3|) . 

a £ = a + O(e) in H s , Vs > 0, 

then the above estimate can be improved to: 

IK - a £ \\ LOO{[ ^ T] . HS) = 0(e), Vs > 0. 

We have reduced the study of the asymptotic behavior of u £ to the un- 
derstanding of a £ . To complete the proof of Proposition 11.21 we resume 
the framework of Section [21 With x(t,y) given by the Hamilton flow (|2.1jl . 
introduce the Jacobi determinant 

J t {y) = deW y x(t,y). 

Denote 

A £ (t,y) :=a £ (t,x(t,y)) s/My)- 
We see that so long as y i— > x(t, y) defines a global diffeomorphism (which is 
guaranteed for t £ [0, T] by construction), ()3.3|) is equivalent to: 

d t A e = -iE K - 1 f(j t (y)- 1 \A B f)A e ; A £ (0, y) = a (y). 

This ordinary differential equation along the rays of geometrical optics can 
be solved explicitly: since / is real- valued, we see that d t \A £ \ 2 = 0, hence 

A £ (t,y) = a (y)exp {-ie^ 1 f (j s {y)- x \a {y)f) ds^j . 

Back to the function a £ , Proposition 11.21 follows, with: 

( 3.4) >M(v(*.*)) 

G(i,z) = - / /(j s (y(t,x))- 1 |a (y(t,x))| 2 )d S . 
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One may wonder if this approach could be extended to some values k < 1. 
Seek a solution of (|3.1|) . To have a simple ansatz as in Proposition 11.21 we 
would like to remove the Laplacian in the limit e — > in (|3.1|) . and obtain 
the analogue of Corollary 13.21 Following the same lines as above, we find: 

a £ (t,x) = a(t,x)e i£K ~ lG{t ' x \ 

which is exactly the first formula in Proposition 11.21 Now recall that in 
Proposition 13.11 we prove that a £ is bounded in H s , uniformly for e e]0, 1]; 
this property is used to approximate a £ by a £ . But when k < 1, a £ is no 
longer uniformly bounded in H s , because what was a phase modulation for 
k > 1 is now a rapid oscillation. 

This remark in a particular case reveals a much more general phenomenon. 
When studying geometric optics in a supercritical regime (when k < 1 in 
the present context), distinguishing phase and amplitude becomes a much 
more delicate issue. Suppose for instance that we seek u £ = a £ e l ^^ £ , where 
a £ and (ft £ have asymptotic expansions as e —* 0. All the terms for <p £ which 
are not o{e) are relevant, since (fi £ is divided by e. To determine these terms, 
it is not sufficient to determine the leading order amplitude lim a £ in general: 
because of supercritical interactions, initial perturbations of the amplitude 
may develop non- negligible phase terms. To illustrate this discussion in the 
above case, let k < 1 and a £ solve ()3.3|) with initial data 

Ou_ = a + e 7 ai, 7 > 0, 

where 00 and a\ are smooth and independent of e. Integrating (|3.3|) . we find 

a e (t,z) = 1 (oo (y(t,x))+sM (y(t,x)))e i£K - lG ^\ 
y/Jt[y(t,x)) 

where 

G £ (t, x) = -J t f (j s (y(t, x))- 1 \a (y(t, x)) + e> ai (y(t, x))| 2 ) ds. 

We have the identity G £ = G + e^Gi +e 2 ' y G2, for the same G as before, and 
G\,G2 independent of e, but depending on 01. We infer: 

a £ {t, x) ; 1 a x)) e i£K " 1 ( G + £7G i+ e27G 2)(< i x) > 

In particular, if K + 7 < 1, we see that G\ has to be incorporated to describe 
the leading order behavior of a £ . Since the three requirements k < 1, 7 > 
and k + 7 < 1 can be met, we see that the small initial perturbation e 7 ai 
may produce relevant phase perturbations. This example explains why in 
Assumption 13 we require the asymptotic behavior of the initial amplitude 
up to order o(e) and not only o(l) (take k = and 7 = 1). We refer to ^2] 
for a general explanation of this phenomenon, and to the next three sections 
as far as Schrodinger equations are concerned. 

4. A HYPERBOLIC POINT OF VIEW 
In this section, we study (|1.1|) in the C8. SC K = 0, with no potential: 



(4.1) ied t u £ + £ -&u £ = f (\u £ \ 2 ) u £ ■ ut =0 = a £ (x)e 
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We recall the method introduced by E. Grenier [22], which is valid for smooth 
nonlinearities which are defocusing and cubic at the origin. Throughout this 
section, we assume the following: 

Assumption 4 (Study of (|4.1|l). We have f > 0. In addition, <p G H°°, 
and there exists ao,a\ E H°° such that 

o-o = o-o + eai + o(e) in H s , Vs > 0. 

Note that this assumption is closely akin to Assumption |2J nevertheless, 
we do not make any assumption on the momenta of ao and a%, and the initial 
phase is bounded. We will see in Section |S] how to weaken this assumption. 

4.1. Grenier's approach. The principle is somehow to perform the usual 
WKB analysis "the other way round". First, write the exact solution as 

(4.2) u £ {t,x) =a £ (t,x)e i * £(i ' x)/e , 

where $ £ is real- valued. Then show that the "amplitude" a E and the "phase" 
<£ £ have asymptotic expansions as e — > 0: 

a £ ~ a + ea (1) + e 2 a (2) + . . . ; $ £ ~ <j> + e0 (1) + e 2 ^ (2) + . . . 

Introducing two unknown functions to solve one equation yields a degree of 
freedom. The historical approach |26| Chap. Ill] consisted in writing 

^ + i|V^| 2 + /(K| 2 )=, 2 ^ ; n= = ^o, 

d t a £ + V$ £ • Va e + ^a £ A<5> £ = ; a £ | t=0 = ag. 

Of course, this choice is not adapted when the amplitude a £ vanishes (see 
|21|). so it must be left out when Oq E L 2 (M n ) in general. The approach 
introduced by E. Grenier consists in imposing: 

d t <S> £ + l\V$ £ \ 2 + f(\a £ \ 2 ) =0 ; $ £ | =^ , 

(4-3) 2 1 

a t a e + V$ e • Va £ + -a e A$ e = i-Aa £ ; a e |. =Oq. 

2 2 t— u 

Before recalling the results of [22], observe that if a £ and <3? £ are bounded in 
some sufficiently small Sobolev spaces uniformly in e, passing to the limit 
formally in (|4.3|) yields: 

dt4> + ~ |V<£| 2 + / (|a| 2 ) =0 ; 4 =0 = ^o, 
(4.4) 2 x 

5 4 a + V0 • Va + -aA0 = ; a| t=0 = a . 

We see that when the nonlinearity is exactly cubic (f(y) = y), (p,v) '■= 
(|a| 2 , V</>) solves the compressible, isentropic Euler equation 

9 4 u + w • + Vp = ; v\ =V(f , 

(4.5) 

+ V • (pv) = ; p\ t=0 = \a \ . 

From this point of view, the formulation (|4.4|) is closely akin to the change 
of unknown function p — ► ,/p introduced in |2H] (see also ^J) to study (|4.5j) 
when the initial density is compactly supported, a situation more or less 
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similar to the present one. Note however that here, a is complex-valued in 
general. 

Introducing the "velocity" v £ = V$ £ , 1)4. 3 j) yields 



t=o 



l*=o ~ °o • 



d t v £ + v £ -Vv £ + 2f (\a £ \ 2 )Re(a £ Va £ ) = ; u £ 

(4-6) i e 

d t a £ + v £ ■ Va e + -a £ V ■ v £ = i^a £ ; a £ 

Separate real and imaginary parts of a £ , a £ = a\ + %a\ . Then we have 

n 

(4.7) 8 t u £ + Mu^djtf = e -Lu £ , 



3=1 



with u e 



and A(u,Z)=J2 A i( U ^ 

3=1 



-A ... 
A ... 

nxn 









2/'ai C 2fa 2 £ t; • 



where /' stands for /'(|ai | 2 + 1 «2 1 2 ) - The matrix A(u,£) can be symmetrized 
by 



(4.8) 



S 



h o 
o -Li 



which is symmetric and positive since /' > 0. For an integer s > 2 + n/2, 
we bound (Sd% u £ , <9"u £ ) where a is a multi index of length < s, and (•, •) is 
the usual I? scalar product. We have 

j t (Sd«u £ , «9> £ ) = (d t Sd«u £ ,d«u £ ) + 2 (Sd t d%u £ , d«M £ ) 
since S is symmetric. For the first term, we consider the lower n x n block: 



-pM/W + Ki 2 )) 



(d t Sd«u £ ,d«u £ )< 
So long as ||u £ ||£oo < 2||g[q||£°°, we have: 

f'H\ 2 + \a £ 2 \ 2 )>mf{f'(y); < y < 4 sup ||a§||£,» ) = <5 n > , 

I 0<£<1 J 

where 5 n is now fixed, since /' is continuous with /' > 0. We infer, 



-ia,(/'(|af| 2 + |a|| 2 )) 



< C||u e || H . 



L°° 



where we used Sobolev embeddings and (|4.7jl . For the second term we use 



(Sd t d«u £ ,d«u £ ) = - (SL(«9> £ ),d> £ ) - [SdS[J2 A ^ £ )dj^),d> £ 

3=1 
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We notice that SL is a skew-symmetric second order operator, so the first 
term is zero. For the second term, use the symmetry of SAj(u £ ) and usual 
estimates on commutators to get finally: 

j t £ (Sd^,d^) < c(kim £ (S8£v?,asvf) , 

\ct\<s M<s 

for s > 2 + d/2, Gronwall lemma along with a continuity argument yield: 

Proposition 4.1 (|221> Th. 1.1). Let Assumption^be satisfied. Let s > 2 + 
n/2. There exist T s > independent of e G]0, 1] and u e = a e e l * e / £ solution 
to 1)4 .lj> on [0,T S ]. Moreover, a £ and <F are bounded in L°°([0, T 8 \\ H s ), 
uniformly in e g]0, 1]. 

The solution to Q4.3JI formally converges to the solution of (|4.4jl . Under 
Assumption 01 ()4.4() has a unique solution (a, <^>) 6 L°°([0, T*]; H m ) 2 for any 
m > for some T* > independent of m (see e.g. (21I1E1)- We infer: 

Proposition 4.2. Let s E N. T/ien T s > T#, and there exists C s independent 
of e such that for every < t < T* , 

(4.9) ||a £ (t) - a(*)|| H « < C s e ; ||$ £ (t) - 0(t)|| ff . < C s et. 

Proof. We keep the same notations as above, 1)4.7(1 . Denote by v the analog 
of u £ corresponding to (a, (ft). We have 

d t (u £ - v) + M" £ )dj (u £ - v) + £ (A,(u £ ) - ^(v)) iV = |Lu £ . 

J=l 3=1 

Keeping the symmetrizer S corresponding to u £ , we can do similar computa- 
tions to the previous ones. Note that we know that u £ and v are bounded in 
L~([0,min(T S) T*)];i2" s ). Denote w e = u e -v. Writing Lu £ = Lw e +Lv, the 
term Lw £ disappears from the energy estimate, and we get, for s > 2 + n/2: 

j f £ (S9> £ ,9> £ ) <C(||u £ ||^,||v||^ +2 ) (^> e ,9> £ ) 

|a|<s M<s 

+ £\\v\\ H s+2\\w £ (t)\\ H s ■ 

Gronwall lemma and a continuity argument show that w e (hence u e ) is 
defined on [0,T*]. By Assumption^! ||wf t=0 ||.H-s = 0(s), and we get: 

ll w£ Hi°°([0,T«];if s ) = ^( £ )- 

The estimate for the phase (and not only its gradient) then follows from the 
above estimate and the integration in time of (|4.3|) - 1)4.4(1 . □ 

Proposition 14.21 yields an approximation of u £ for small times only: 

\\u £ (t)-a(t)e^ it)/£ \\ L2 = ||a E (t)e i<I>e(i)/e - a(t)e^ (i)/e 

= 0(\\a £ (t)-a(t)\\ L2 + e ^(t)/e_ e m)/e \\ a{t )\\ L2 

= 0(e) + 0(t). 

For times of order 0(1), the corrector a\ must be taken into account: 
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Proposition 4.3. Let Assumption^ be satisfied. Define (a^,^ 1 )) by 
d t ^ 1] + • Vc/> (1) + 2 Re (aa (1) ) /' (|a| 2 ) = 0, 



d t a {1) + V0 • Va (1) + Vc/> (1) • Va + -a (1) Ac/> + -aA0 (1) = ^Aa, 



^(i)| =Q . (i)| 
V \t=o u ' " \t=o 



at. 



Then a^\(p^ G L°°([0, T*]; H s ) for every s>0, and 

||a £ -a-ea (1) || L oo ([OiTsfc ].^ s) + ||$ e -0-e(/) (1) || L o O([OiTsfc ] ;Hs) < C s e 2 , Vs > . 

The proof is a straightforward consequence of the above analysis, and is 
given in [22]. Despite the notations, it seems unadapted to consider as 
being part of the phase. Indeed, we infer from Proposition 14.31 that 



u £ -ae^ (1) e^ /e 



= C(e) 

L°°([0,T«];L 2 nL°°) 



Relating this information to the WKB methods presented in the introduc- 
tion, we would have: 

a = ae t(p . 

Since 4>^ depends on a\ while a does not, we retrieve the fact that in super- 
critical regimes, the leading order amplitude in WKB methods depends on 
the initial first corrector a\. 

Remark 4.4. The term e*^ (1) does not appear in the Wigner measure of 
ae 1 ^ 1 ' ] e l ^l e . Thus, from the point of view of Wigner measures, the asymp- 
totic behavior of the exact solution is described by the Euler-type system 

Remark 4.5 (Introducing an isotropic harmonic potential). The above method 
makes it possible to consider the semi-classical of when V(t, x) = \\x\ 2 
is an isotropic harmonic potential, and Assumption |I] is satisfied. Let 

U £ (t,x) = Kr-ne T +t I ^ u e I arctant, ; . 

Then U £ solves: 

ied t w + £ ^au £ = -L./ ((i + t 2 f 2 \U £ \ 2 ) u £ , 

U £ (0,x) =a £ (x)e i ^ x)/£ . 

We can then proceed as above. The only difference is the presence of time 
in the nonlinearity, which changes very little the analysis. 

Remark 4.6 (Momenta). If in Assumption 01 we replace H s with 

S s = ff s nf(ff s ) = {»eL 2 ; (1- A) k / 2 (x) s - k w £ L 2 , < k < s} , 

then the above analysis can be repeated in X s . The main difference is due 
to the commutations of the powers of x with the differential operators; it is 
easy to check that they introduce semilinear terms, which can be treated as 
source terms when applying Gronwall lemma. 
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4.2. Remarks about some conserved quantities. Consider the case of 
the cubic, defocusing Schrodinger equation: f(y) = y. Recall three impor- 
tant evolution laws for 

Mass: ^\\u £ (t)\\ L 2 =0. 
Energy: - (\\eV x u £ \\ 2 L2 + \\u £ \\ A Li ) =0. 
Momentum: — Im / u £ (t, x)eX7 x u £ (t, x)dx = . 



Pseudo-conformal law: {\\J e (t)u £ \\ 2 L 2 + ^ll^lll-O = *( 2 _ n )ll ne |ll 4 ' 



where J £ (t) = x + iefSJ x . These evolutions are deduced from the usual ones 
(e = 1, see e.g. ^QJ El]) via the scaling i^{t,x) = u{et,ex). Using (|4.2|) and 
passing to the limit formally in the above formulae yields: 



j t \\a(t)\\ L i = . 



dt 



d 
dt 



(|a(t,x)| 2 |V</>(t,x)| 2 + \a(t,x)\ 4 ) dx = . 
d f 

— / \a(t,x)\ 2 V<p{t,x)dx = 0. 

(x - tVcf)(t, x)) a(t, x)\ 2 + t 2 \a(t, x)\ 4 ^jdx = 

= (2-n)t J \a(t,x)\ 4 dx. 

Note that we also have the conservation ([S]): 
d f 

— Re / W(t,x)J £ (t)u £ {t,x)dx = 0, 



which yields: 

d 
di 



((x - tV(p(t, x)) \a(t, x)\ 2 ) dx = . 



All these expressions involve only (|a| 2 ,V<^>), that is, the solution of (|4.5|) , 
We thus retrieve formally some evolution laws for the Euler equation. 



5. Introducing subquadratic potential and initial phase 

In this section, we prove Theorem ll.3l First, we point out that the unique- 
ness for u £ in C([0, T*]; H s ) is straightforward for s > n/2. We thus have to 
prove that there exists such a solution, and that it is smooth. 

As suggested by the statements of Theorem II. M| the idea consists in re- 
suming Grenier's method, and in writing the phase & £ as 
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We take 4> £ as a new unknown function. Recall that the system Q4.3|) reads, 
with the present notations: 



d t § £ + -\V§ £ \ 2 + V + f{\a £ \ 2 ) =0 ; $ 

d t a £ + V$ £ • Va £ + -a £ A$ £ = i-Aa £ ; a £ 
2 2' 

This system becomes, in terms of <f> £ , and given Q1.5|) : 



lt=o 



(5.1) 



d t <f + - |V<f | 2 + Vc/> cik • Vc/> £ + / (|a £ | 2 ) 
d t a £ + V0 £ • Va £ + V0 cik • Va £ + ^a £ A0 £ + ^a £ A0 cik 



0, 



z|Aa £ , 



it=o 







it=o 



a, 



o • 



Like in Section l4,l( we work with v £ = V<j) £ instead of 4> £ , to begin with. 
The new terms are the factors where </> e ik is present. The point is to check 
that they are semilinear perturbations, which can be treated as source terms 
in view of Gronwall lemma. Again, separate real and imaginary parts of a £ , 
a £ = a\+ ia 6 ,, and introduce: 



u 



L 



-AO ... 
A ... 

nxn 



n 

and A(u,f) =J2 A j( u )ti 
i=i 



v-C 

v£ f*c 



21 <f 
2 S 
«2 t< 



where /' stands for /'(|ai| 2 + | a,2 1 2 ) - Instead of (|4.7I) . we now have a system 
of the form 

n n 

(5.2) ^u £ + ^^(u £ )^u £ + ^i?,(V ( /> eik )5 j u £ + M (V 2 </> eik ) u £ = £ -Lu £ , 



where the matrices Bj depend linearly on their argument, and the matrix 
M is smooth, locally bounded. The quasilinear part of (|5.2j) is the same 
as in Section 14.11 and involves the matrices Aj. In particular, we keep the 
same symmetrizer S given by (|4.8jl . The matrices Bj have a semilinear 
contribution, as we see below. The term corresponding to the matrix M 
can obviously be considered as a source term, since ^ e ik is subquadratic. 

Let s be an integer, s > 2 + n/2, and let a be a multi index of length < s. 
We have: 

j t (Sd> £ , d> £ ) = (d t Sd«u £ , <9> £ ) + 2 (S^> £ , 9> £ ) 
since 5 is symmetric. For the first term, we consider the lower rax n block: 



(5.3) (d t Sd%u s ,d%u e ) < 



jjd t (f'(\alf + \alf)) 



(Sd> £ ,9> £ 
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We consider times not larger than T given by Lemma ll. 11 so that the func- 
tion (/^ik remains smooth and subquadratic. So long as ||u e ||ioo < 2||ao||z,°°, 
we have: 

/'(l«i| 2 + l« e 2| 2 ) >inf{/'(y) ; < y < 4 sup ||a§||ioo ) = S n > . 

I 0<£<1 J 



We infer, 



(5.4) 



/ 



<C(||u B || H . + ||a;u 6 || H .-i; 



l CO 



for some locally bounded map C(-). We used Sobolev embeddings, 1)5.2(1 
and Lemma ll.lt the terms Bj are sublinear in x, hence the norm ||xu £ ||jy s -i 
which we did not consider in Section [4.11 We emphasize that this estimate 
explains why we assume s > 2 + n/2, and not only s > 1 + re/2: we control 
dtu £ in L°° using (|5.2[) . so we need to estimate Lu £ in L°°. For all the other 
terms, s > l+ra/2 would suffice. This also explains why we wrote ||xu £ ||jy s -i 
and not ||xu e II h-s. For the second term we use 



n 

i=i 

n 

-(^(E^( V ^)^ u£ )' a > £ ) " (^(M(vVcik)u £ ),9> 



The first two terms of the right hand side are handled in the same way as 
in Section [4.11 the first one is zero, and the second can be estimated by: 

n 

(5.5) (S^( J> i (u e )a / u e ),^u e ) < C(\\u £ \\ Hs ) (Sd«u £ ,d«u £ ), 

j=l |a|<s 

where we keep the convention that C(-) is a locally bounded map. Let us 
briefly explain this quasilinear estimate. First, write 

(S3" (^(u'jaX) , d> £ ) = (SMn'fadZvf, d> £ ) 

+ (5 (^(^(u £ )^u £ ) - A,-(u e )^u 6 ) , 9> £ ) . 

By symmetry of SAj(u £ ), 

(SAiWdjdZxf&vf) = - (^(^(u £ ))3> £ ,9> £ ) 

-(SA 3 (u £ )d 3 d«u £ ,d«u £ ). 

We infer: 

|(SA> £ )^> e ,d> £ )| < WdjiSAj^))^ \\d«u £ \\ 2 L2 



< C(||u E || roo ')||V 1B U 6 ' 



| £jOO I || V x u || £00 

The usual estimates on commutators (see e.g. j^S]) lead to 

|(5 (A,(u £ )<9,u £ ) - A^u^-^u 6 ) ,^u 6 )| < C(||u £ || Hs ) ||u e ||^ , 
and (|5.5|1 follows, since we consider times where S _1 is bounded. 
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For the third term of (Sdtd"u £ , d^u £ ), write: 

+ J S [d^BjiV^d^^d^dx. 
For the first term of the right hand side, an integration by parts yields: 
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(5.6) 



SBj (yMd&vfdZvFdx 



< WdjiSBjiV^W^wu 



er || 2 



(x) a £ Va £ 



< C(||u e || L oo)Q|u e || LOO + ||xu e || LOO + ||Vu e 



u 



|u £ ||^ 



e II 2 



where we have used Lemma ll.ll Again from Lemma the commutator 

[d^BjiV^dj] 

is a differential operator of degree < s, with bounded coefficients. We infer: 



5^(^^(V<Aeik)^u £ j,5> e J < C(||u £ ||^ + ||xu £ || Hs -i)||u £ ||^ 
i=i 



We have obviously 



S^(M(V^eikK),d> £ ) < C(\\u £ \\ L ~)\\u £ \\ 2 Hs . 



This yields: 

(5.7) j t {S8gxf,8£u<) < C(\\u £ \\ Ha + \\xu?\\ H ^) \\u £ f Hs , 

where the map C(-) is locally bounded. We now have to bound xu £ in iiP -1 
to close our family of estimates: we consider 

'sdf (x k u £ ), d? (x k u s )) ; 1 < k < n, \(3\ < s - 1. 



dt 

We can proceed as above, replacing u £ with x k u £ : x k u £ solves almost the 
same equation as u £ , and we must control some commutators. 

n n 

d t (x k u £ ) + J2M" £ )dj(xkU £ ) +^B J (V4 ik )a j (^u e )+ 

3=1 3=1 

+M (V 2 0eik) x k u £ = ^L{x k u £ ) + A k (u £ )u £ + B k (V<p cik )u £ + |[x fcj L]u £ . 

The term A k (u £ )u £ is harmless. The term i?fc(V0 o ik)u £ is controlled by 
(x) u £ , since e ik is subquadratic: this is a (semi)linear perturbation. Finally, 

/ 2d k ... 
[x k ,L] = -28 k ... 

V o o o nxn 

Now we only have to notice that estimating x k u £ does not involve extra 
regularity or extra momenta. In the above computations, the first time 
we needed to consider momenta was for (|5.4|) : we need exactly the same 
estimate now, since it is due to the symmetrizer, which remains the same. 
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The same remark is valid for (|5.6[) . For j3 a multi index of length < s — 1, 
we find: 



(5.8) M -M £ ).^H< 

< C7 (||u £ ||^ + Wxu'WB-i) {\\u £ \\ 2 H s + II^H^) . 

The term ||u e ||^ 3 is due to the commutator [xfc,L]: 

'Sd% ([x fc ,L]u £ ),df (x k u e )) = fe([x fc ,L]u £ ),d£(x,u £ 



< ||u e ||#s ||xu e ||^ s -i. 

Summing over the inequalities (|5.7|) and (|5.8() yields a closed set of estimates, 
from which we infer the analogue of Proposition 14.11 note that the time T s 
is not larger than T by construction, and may be strictly smaller than T, 
due to possible shocks for (|5.2|) . We also mention the fact that the above 
analysis gives v 6 = V<j) 6 £ C([0, T s ]; H s ), with xV(j) e £ C([0,T s ]; H 8 ' 1 ): 
back to <GZH), this shows that 4> £ G C([0, T*]; H 3 ' 1 ), but that we cannot 
claim that x<p £ £ C([0, T*]; iP -1 ). 

Passing to the limit e — ► in Q5.1|) . it is natural to introduce the system: 

£W + \ I V0| 2 + V0 cik • + / 0| 2 ) = o, 

( 5 - 9 ) d t a + V0 • Vo + V4 ik • Va + -aA<j) + ^aAc/> eik = 0, 

^t=o = ! a \t=o = a o- 
The above analysis shows that this system has a unique solution in H s 
with the first momentum in iJ s_1 , locally in time for t G [0,T*], for some 
T* e]0,T]; T* is independent of s, from the usual continuation principle, 
explained for instance in |28[ Section 2.2]. We easily obtain the analogue 
of Proposition 14.21 mimicking the above computations, we can estimate the 
error (a e — a, eft 6 — <fi) in H s , by first estimating (a £ — a, V(f) e — Vcf)) and its 
first momentum in H s . We deduce that T s > T*, and (|1.8|) follows. 

The end of Theorem II. HI can then be proved as in : from the above 
analysis, the functions 

a 6 - a V(j) e - V4> 

£ £ 

and their first momentum, are bounded in H s for every s > 0. A subse- 
quence converges to the linearization of (|5.2|) . yielding a pair 
By uniqueness for the limit system, the whole sequence is convergent, and 
the analogue of Proposition 14.31 follows. This completes the proof of Theo- 
rem ESI 

6. Extension to the case < k < 1 

When < k < 1, we propose an analysis which can be considered as a 
generalization of the study led in Section Throughout this paragraph, 
we suppose that Assumptions lll il are satisfied. Again, we write the exact 
solution as 

u e = a¥^ £ , with $ £ = 4ik + <f- 
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The unknown function is the pair (a £ , (f) £ ). We have two unknown functions 
to solve a single equation, We can choose how to balance the terms: 

we resume the approach followed when k = 0. Note that this approach 
would also be efficient for the case n > 1, with the serious drawback that 
we still assume /' > 0, an assumption proven to be unnecessary when k > 1 
(see Section EJ). We impose: 

dtf + \ |V<f | 2 + V0 eik • V0 £ + e K / (|a £ | 2 ) = 0, 

( 6 - 1 ) d t a £ + V<f ■ Va £ + V0 cik • Va £ + ^a £ A0 £ + ^a £ A</> cik = i|Aa £ , 

<f| t=0 = ; a £ \ t=0 = a £ . 

This is the same system as with only / replaced by e K /. Mimick- 

ing the analysis of Section 03 we work with the unknown u £ given by the 
same definition: it solves the system 1)5 .2 J) , where only the matrices Aj have 
changed, and now depend on e. The symmetrizer is the same as before, 
with /' replaced by e K /': the matrix S = S e is not bounded as e — > 0, but 
its inverse is. We see that 1)5. 3 Jl and 1)5.4)1 still hold, independent of k. We 
claim that inequalities similar to (|5.7jl and (|5.8jl hold: 

| (5 £ 9> £ ,0> £ )<C(||u £ || HS + llxu-ll^-x) £ (S £ d> £ ,d> £ ), 
| (Vdf(* fe u £ ),df (x fc u £ )) < C(||u £ || HS + \\xn<\\ H .-i)( E (^> £ ,<9> £ ) 

l<j<n 

where the map C(-) is locally bounded, independent of e g]0, 1]. The fact 
that such estimates remain valid, with this dependence upon e, stems essen- 
tially from the following reasons: 

• The matrices S £ and Bj are diagonal. 

• The matrix M is block diagonal: the blocks correspond to the pres- 
ence/absence of e in S £ . 

• The matrices S £ Aj are independent of e e]0, 1]. 

• The inverse of S £ is uniformly bounded on compacts, as e —* 0. 

A continuity argument and Gronwall lemma then imply the analogue of 
Proposition 14.11 u £ exists locally in time, with .£P-norm uniformly bounded 
as e — > 0. Note that since £ | t _ o = 0, we have: 

(5 £ 5> £ ,9> £ )| i=0 = O(l), 

and we infer more precisely: 

IKIU°°([0,T,];H*) = ; \\V4> e \\L°°{[0,T*\;H°) = ° O 5 ^) • 

It seems natural to change unknown functions, and work with <p £ = e~ K <j) £ 
instead of <j) £ . With this, we somehow correct the shift in the cascade of 
equations caused by the factor e K in front of the nonlinearity. Then 1)6.1)1 
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becomes: 



V0 £ 



+ V0 eik - V0 £ + /(|a £ | 2 ) =0 



( 6 - 2 ) <V £ + e K V</> £ • Va £ + V^ik ■ Va £ + y a £ A</> £ + ^Afe = i|Aa £ , 



\t=Q 







it=o 



a6 



The pair (<^> £ ,a £ ) is bounded in C([0,T]; H s ). Therefore, a subsequence is 
convergent, and the limit is given by: 



dt<f> + V<t>cik- V0 + /(|a| 2 ) =0 ; 
d t a + V(?!> e ik • Va + -aA0 cik = ; a| 



li=o 



t=o 



0, 



a 



We see that a solves (|l.f)[) : is given by an ordinary differential equation 
along the rays associated to ^ e ik> with a source term showing nonlinear 
effect: / (|a| 2 ). By uniqueness, the whole sequence is convergent. Roughly 
speaking, we see that if 



w 



V 



then Gronwall lemma yields: 

(S £ d£w £ , d£w £ ) < C (e + e K ) < 2Ce K . 

We infer: 

Proposition 6.1. Let s > 2 + n/2. Then (|6.1|) has a unique solution 
(a £ ,0 £ ) e C([0,T];iT s ) 2 ; smc/j i/wrf x k a e , x k dj(j) £ £ C([0,T]; H s ), for every 
^ ^ j,k < n (T is given by Lemma \l.l\) . Moreover, there exists C s indepen- 
dent of e such that for every < t < T , 



(6.3) 



: (t)-e K <p\\ H s <C s e z H 



\\a £ (t)-a(t)\\ H s <C s e K ; 
where a is given by (jl.fi j) . 

Three cases must be distinguished: 

• Ifl/2<K<1, then we can infer the analogue of (|1.9|) . 

• If k = 1/2, then we can infer the analogue of (|1.8|) (but not yet of 

(EHD). 

• IfO < k < 1/2, then we must pursue the analysis, and compute a 
corrector of order e 2k . 

We shall not go further into detailed computations, but instead, discuss 
the whole analysis in a rather loose fashion. However, we note that all the 
ingredients have been given for a complete justification. 

Let N = [1/k], where [r] is the largest integer not larger than r > 0. We 
construct a^ 1 ', . . . , and f^ 1 ), . . . , 0^ such that: 



+ 



L°°([0,T];H s ) 



+ 



L°°([0,T];H s ) 
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But since N + 1 > we have: 



,2^(1) _ ... _ £ N«ftN- 



The analogue of ()1.9|) follows: 



i) 



L°°([0,T];i/ s ) 



/ £ (iV+lK\ =o(£) 



where 



ae 



app 



L°°([0,T];L 2 nL°°) 



0(1), 



rl-K 



+ 



gl) 

-1-2/c 



+ 



+ 



£V-F) 
-l-JWc ' 



Remark 6.2. In the case « = 1, iV = 1, and the above analysis shows that 
one phase shift factor appears: we retrieve the function G of Proposition E21 
(under the unnecessary assumption /' > 0). If k > 1, then N = 0, and we 
see that ae l ^ cik ^ e is a good approximation for u e . 
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